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Abstract 



We reexamine the recent instanton motivated studies of Alford, Rajagopal 
and Wilczek, and Berges and Rajagopal in the framework of the standard 
SU(2) Nambu-Jona-Lasinio model. The chiral phase diagram is calculated in 
the temperature-density plane, and the pressure is evaluated as the function 
of the quark density. Obtaining simple approximate relations describing the 
T-fj, and T-pp phase transition lines we find that the results of the instanton 
based model and that of the NJL model are identical. The diquark transition 
line is also given. 
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I. INTRODUCTION. 



Recent studies by several authors using an effective 4-fermionic interaction between 
quarks []T]-§J or direct instanton approach || have rekindled interest in the two flavor QCD 
phase transitions. In particular, Alford, Rajagopal and Wilczek |IJ have studied the pres- 
sure density and gap parameter using a fermionic Lagrangian with an instanton motivated 
four-point interaction. At zero temperature these authors found negative pressure for a cer- 
tain range of the Fermi momentum, p-p and showed the solutions of the gap equation as a 
function of pp. Berges and Rajagopal [0] extended this work and have calculated the phase 
diagram of a strongly interacting matter as a function of temperature and baryon number 
density in the same model. 

The question that we raise and examine in this paper is whether or not these results 
are fundamentally different from those obtained via the standard well-known Nambu-Jona- 
Lasinio (NJL) model [§H9j, at least with regard to chiral symmetry breaking. Quantities 
such as the gap parameter, pressure density and other thermodynamical quantities have 
been extensively studied in this model over the last decade p|,[7|, p!0|JTT[| , and even to a level of 



sophistication that goes beyond the standard mean field treatments [10,11|. However, the 
results have usually been presented as a function of the chemical potential, and not as of the 
Fermi momentum pp or density, as the authors of have done, and hence the connection 
between their results and those of the NJL model are not obvious. Thus, in order to make 
a systematic comparison, we have to reevaluate the gap, pressure and phase diagram in 
these variables. In section II we derive a simple analytical approximate expression for the 
phase boundary that is independent of the parameters of the NJL model. In section III, we 
examine the pressure as a function of the density and calculate the complete chiral phase 
diagram numerically from Maxwell constructions. We compare our results with that of Ref. 
0. In Section IV, we write down the form of the gap equation for a superconducting diquark 
transition, and thus the line of points in the T-/J, plane. This is evaluated numerically. We 
summarize and conclude in Section V. 



II. PHASE BOUNDARY CURVE - AN ANALYTIC EXPRESSION. 

We commence by deriving an approximate analytic expression for the phase boundary 
curve. Our starting point is the gap equation for the dynamically generated up and down 
quark mass m that is derived from the SU(2) chirally symmetric Lagrangian, 

£njl = $W + G[(^) 2 + ($75 r^P) 2 } (1) 
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with G a dimensionful coupling and ip quark spinors for u and d quarks. At non-zero 
temperature and chemical potential, the mean-field self-energy or gap equation reads || 

S* = m = AGN c N f m J-^^[l- - r (?,!*)] (2) 

with the Fermi distribution functions = \/\e^ Ep±,x > + 1] = f ± (p,fj,). The gap equation is 
easily seen to minimize the thermodynamical potential p|,|Tl|, 

fl(m) = g - 7 / -irj-j^L log[l + e~^} [1 + e~^\ (3) 

where E 2 = p 2 + m 2 , (3 = 1/T and 7 = 2N c Nf is the degeneracy factor. The condensate 
is related to m via m = —2G(ipip). The three momentum integrals are understood to 
be regulated by a cutoff A, and a standard set of parameters, A = 0.65 GeV and G = 
5.01 GeV -2 are used to fix the values of f n = 93 MeV and the condensate density per flavor, 
(uu) = (dd) = (-250 MeV) 3 . 

Let us examine first the T = limit of Eq.(@). The gap equation for the non-trivial 
solution is 

71-2 -a p2 
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in terms of the Fermi momentum pp. p-p is a decreasing function of the constituent mass 
taking its maximum value, p c at the chiral phase transition, m — > 00 



Pc = A A 



7T 2 



(5) 
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Using the numerical values of G and A given above and fixing N c = 3, Nf = 2 leads to 
the numerical value p c = 0.307 GeV. We use this relation to eliminate A 2 in favor of p 2 . 
Now, in order to determine the phase transition line, we start with the gap equation (0) and 
divide m 7^ out. The critical values of the temperature and chemical potential that lie 
on the phase transition boundary are determined through the vanishing chiral condensate 
and hence m — > 0. Performing the usual substitution u = /3(p ± fi) as appropriate, and 
evaluating the integral on the right hand side, one finds 



1 Here we assume a second order transition where the smooth vanishing of the order parameter 
signals the transition. However, for small temperatures the phase transition is first order with a 
jump in the order parameter. Nevertheless, the approximation made here qualitatively gives the 
proper answer (cf. Fig. 2.). 
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This equation defines the chiral phase transition curve in the T-/i plane, and in this form 
contains no explicit model dependence on G and A. We would however, prefer to display the 
phase transition line as a function of temperature and quark density. To do so, we express 
the density through the Fermi momentum, p-p at zero temperature, 



2AUV 
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At finite temperature, the quark density n is given as 

d 3 p 



n 



2N c Nf 



(2k)* 
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(7) 



(8) 



The cube root of the density defines n 3 as n 1//3 
to be 



cn 3 , where c = 2/n 2 , and n 3 is evaluated 



n 3 



3 sinh 



dp 



P 



cosh (/i/T) + cosh (E p /T) 



1/3 



(9) 



Using this transformation, we can display all quantities as functions of n 3 (or density) and 
not of the chemical potential, /i any more. One can take a set of T and [i and calculate for 
example m(T, /i) and n 3 (T,fi). The result obtained can be seen in Fig. 1 where for fixed 
T, T = the mass m is plotted as a function of pp (solid line) and of /i (dashed line) for 
comparison. We note that the behavior found is qualitatively the same as that in 0. At 
the chiral phase transition point, m — >• and n 3 — > n 3)C with 



n 



3,c 



(10) 



This expression can be used to eliminate \x from Eq. @ and the chiral phase transition 
curve in the T-n 3 plane is determined from the simple analytical expression, 

4tt 6 
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T 6 -vr 2 ^T 2 + (^ c -^)=0. 



(11) 



The solution of this third order equation in T 2 , defining the critical temperature, T c at 
which the phase transition occurs for a given ratio p = n 3 ^ c /p c can be written as 



(12) 



with 
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1 < p < 2 1 / 6 



(13) 
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The trivial solution to the gap equation is the only solution if the temperature exceeds the 
maximum value of the critical temperature on the transition curve, 

V3pc 



T 



0.169 GeV 
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or if ris ,c/p c > 2 1/6 . One can now recast the solution for n^ )C in terms of T m yielding a 
parameter-free normalized relation, 



, Pc 
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(15) 



We illustrate this relation in a more usual way in Fig. 2, plotting the critical value of the 
temperature as a function of the quark density divided by normal nuclear matter density 
n = 0.17fm" 3 . 



III. PHASE TRANSITION CURVE VIA MAXWELL CONSTRUCTION. 

We return now to the thermodynamical quantities. From Eq.(|3]) for Q(m), it follows that 
the pressure density is given as 

v = -Q = 1 f ^ In \(l + e-^+^Vl + e -PiE P -»)Y\ , f d3 P E 
The energy density is found to be 

«=VlH* + tfW|^ + S- (17) 

In the limit T— >0,/i— >0, one has 

m* 2 N c N f f A „ r~ r , , 

e vac = e T ^ 0l/ ^o = 777 t~ \ dpp ^p 2 + m*\ (18) 

4G 7T Z J 

where m* = m(T = 0,/i = 0), while an analogous calculation for the pressure density p at 
T — > 0, p — ► yields 

Pvac ^vac • (1^) 

p vac and e vac are independent of temperature and chemical potential and their value is (up 
to the sign) the same. For our choice of parameters, e vac = — (407MeV) 4 . Measuring the 
pressure and energy densities relative to their vacuum values, we have 
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Cphys — £ ^vac j (20) 

Pphys = P + e vac . (21) 

In Fig. 3, we plot the pressure density as a function of 713 for a range of temperatures. Note 
that at T = 0, p becomes negative and displays a cusp like structure. This is bought about 
by the fact that the two solutions for the gap equation enter into the pressure density on 
the different arms of the pressure: the rising curve is bought about by the m = solution, 
while the solution that goes down has a value of m 7^ 0. Note that this situation is similar 
to that observed in Ref. |l|]. 

The phase transition curve can now be calculated using a more standard but numerical 
treatment and compared to the approximate (second order) curve shown as the dashed line 
in Fig. 2. The difference between these two curves lies in the fact that the analytic one is 
calculated with A —>■ 00, while the numerical curve has A finite. It is worth noting that 



the truncation of the NJL model to zero modes |L2| shows a similar behavior to the one 
observed in the full model, however, lacking the unstable backbending of the phase curve at 
high temperatures. Plotting the pressure rather as a function of volume allows one to perform 
Maxwell constructions and obtain the full information of the phase diagram, including the 
metastable region. The results of this calculation are show in Fig. 4 as a function of n 1//3 
as the solid curves. Also shown (dashed curves) are the calculations of Ref. [0]. We note 
that the phase transition curve for the chiral transition are qualitatively identical. The main 
difference between these curves lies in the observation that the mixed phase given by Ref. |2j 
starts already at n — 0, similarly to the truncated NJL calculation ]n|. The equivalence of 
these models is perhaps not simply apparent when one examines the Lagrangians. However, 
the thermodynamical potential of Ref. is precisely that of Eq. (|3|) in the absence of the 
diquark condensate. Thus the differences observed in Fig. 4 can only be attributed to the 
use of slightly different parameters, plus the different method of implementing regularization 
- in the NJL model here, a 'hard' cutoff A is employed, and in the approximate expression 
for the phase curve, is also taken as A — > 00, while the authors of [|l],[| use a soft form factor 
F(p) = A 2 /(p 2 + A 2 ) to regulate their momentum integrals. The physics however cannot 
and does not depend on this and the qualitative results remain unchanged. 

Limitations and difficulties of this model in describing thermodynamical quantities are 



well-known [11,13]. We illustrate one problem in showing the energy per quark plotted at 
T = as a function of the density in Fig. 5. This quantity does not posses a minimum 
at normal nuclear matter density as expected, in apposition to recent linear sigma model 



calculations HT5 
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IV. DIQUARK CONDENSATE TRANSITION LINE 



The thermodynamic potential for three colorsQ Q(m) given in Eq.@ corresponds pre- 
cisely to the A = limit of the function Q(m, A) || 

% A) = g + ^ - 2N f J ^{(N c - 2)[E P + Tln(l + e^ E ^)(l + e^ E ^)] 

+ y/tl + A 2 + s^ 2 + A 2 + 2Tln(l + e^U* 2 )^ + e -^V«-+ A2 )}, (22) 

where £± = _E p ± //, and A is the color superconducting condensate in the diquark channel. 
s is a sign function, and s = ±1 for E p > fi or E p < fi. Here the expression from has been 
modified to remove the form factor - instead a 3-D cutoff \p\ < A as is usual in the NJL 
model and which was used in the previous section is to be understood. G\ is a new coupling 
strength mediating a four-fermion interaction that is attractive in the diquark channel. The 
form for Q(m) = Q(m, 0) from Eq.(|3|) thus indicates that the models are the same, and the 
same gap equation for the chiral transition is retained. For the superconducting sector, the 
gap obtained by differentiating Q(0, A) with respect to A is 

-L = 2N F [ ; 1 — tanh^/a + A 2 + ^^_tanh^/e + A 2 }. (23) 

2d F J (27r)3l v /^^ 2^+ v /e + A^ ^ $ 

This expression is a relativistic generalization of the superconducting gap equation for elec- 



tron pairs [16] in which the quasiparticle energies w£ 2 + A 2 and + A 2 relative to the 



Fermi surface are introduced. Note that at T = 0, one recovers the result of [|TJ for \i ^ 0, 
assuming that the form factor of these authors is set to one. The phase transition line in the 
T-fi or T-n 3 planes can now be obtained quite simply. Assuming that the superconducting 
phase transition can only occur in the region where chiral symmetry is restored, we may set 
E = p. Furthermore we assume the transition to be second order, driven by the condition 
A — > 0. Thus the T-\x critical curve satisfies 

r 2 M 
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f A 1 1 f A 1 1 

/ dpp 2 tanh ~/3(p + u) + / dpp 2 tanh -/3(p — u) 

Jo v + U 2 J a v — U 2 



2NfG\ Jo p + n 2 Jf_i p — n 

/•Mil 

+ / dpp 2 tanh -/?(//-/>), (24) 

Jo n — p 2 

which, with obvious changes of variables reduces to 

7T 2 M+a* u 2 1 M- a« a 2 1 

— — = / dat -2fi+^-) tanh -ft + / d^ + 2fi + ^-) tanh -ft. (25) 
iiVjGi Jo q 2 Jo t, 2 



2 for two colors a massless pionic diquark may form. 
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Unlike the case for electron pairing, in which fi ^> uj d, with up the Debye frequency, we 
cannot regard the logarithmic term as being leading, and it is not possible to obtain a simple 
analytic expression for the right hand side of Eq.(B^). We thus solve Eq.(25) numerically 
for the critical line. Clearly this depends on the choice of the strength G\ and is a sensitive 
function thereof. Arbitrarily demanding that T c = 40MeV at fi — 0.4GeV close to the 
values of Ref. |2| sets G\ = 3.10861 GeV -2 , or G\A 2 = 1.31. The resulting curve is indicated 
by the dotted line in Fig. 4. As can be seen, the qualitative behavior of the model of 
Ref. H is confirmed. We found a somewhat lower gap, A being ~ 35MeV at the chiral 
transition point and zero temperature increasing up to ~ 95MeV at and /i = 0.53GeV and 
zero temperature. The general behavior of the gap parameter as a function of the chemical 
potential and temperature is given in Fig. 6. 

Finally, we comment that although the diquark phase transition line was investigated 
here under the expectation that chiral symmetry is restored, this is not necessarily the case: 
in principle, the diquark phase transition line can extend into the region in which chiral 
symmetry is broken, i.e. where (ipip) ^ or m ^ 0. In practice, this turns out to be a 
function of the parameters chosen. If the diquark phase transition line enters into the region 
of chiral symmetry breaking at a temperature larger than the tricritical temperature, the 
dependence of m(T, jj) that enters into Eq.(23) is continuous, and a solution to this equation 
can be found. For our choice of Gi, however, the diquark transition line would enter into the 
region where a first order phase transition takes place. Thus m(T, /z) is discontinous, and no 
physically accessible solution to this equation can be found. The situation is illustrated in 
Fig. 7, in which Eq.(23) is inverted and solved for G\ and the functional dependence, as well 
as the numerical parameter value 3.1 GeV -2 are plotted as a function of temperature for 
several values of fi. One sees from this figure that the line G\ = 3.1GeV~ 2 cannot intersect 
the curve /1 = 0.3GeV for example, which lies in the broken phase, and therefore there is no 
physically attainable solution. Note however that by adjusting the value of the constant G± 
to a somewhat smaller value could admit a solution within this region. 



V. CONCLUSIONS 



In analyzing the chiral phase transition in the NJL model at finite temperature and 
density, we find the same behavior for the chiral and diquark phases as that reported in 
Jl|] and ||, which use an instanton motivated interaction that is also four point in nature. 
That this must occur can be seen directly from the explicit form of the thermodynamical 
potential that is well-known in our case 0,0, and which is obtained from M on setting 



S 



the form factor to one and introducing a 3-D cutoff A. In addition, we are easily able to 
give an approximate analytic form for the chiral phase curve in the T-fi and T-n^ planes 
that is independent of the model parameters. We have examined the extended form of the 
thermodynamic potential that makes provision for a diquark condensate and obtained the 
appropriate critical line in the NJL model. Our qualitative results conform with those of 
0. In addition, we find evidence for the appearance of a diquark condensate also within 
the region where chiral symmetry is not restored, but this is strongly parameter dependent. 
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Figure Captions 



Figure 1: The gap parameter m, shown both as a function of /i (solid curve) and pp (dotted 
line) at T = 0. 

Figure 2: The phase diagram calculated from the approximate analytically form Eq.flTTD 
(dashed line) and as determined numerically (solid lines) as a function of the quark density 
n, scaled by normal nuclear matter density no = 0.17 fm -3 . 

Figure 3: Pressure shown as a function of 77,3 = (2/ir 2 )n (left graph) and as a function 
of the volume, normalized by a characteristic volume. 

Figure 4: Direct comparison of the NJL phase diagram (solid lines) shown as a func- 
tion of n 1 / 3 , with that of (dashed lines). The lines to the right of n 1 ^ 3 ~ 0.19 are the 
superconducting transition lines. 

Figure 5: The energy per quark, shown as a function of n/no. The minimum occurs 
at approximately n ~ 5n and not at n = 3no 

Figure 6: The diquark gap parameter A, shown as a function of /z and T. The con- 
tour values are given in MeV. 

Figure 7: The functional dependence of Gi as calculated from Eq.(23) is plotted as a 
function of the temperature for different values of the chemical potential. The curves, taken 
from the uppermost one, correspond to the values of /i = 0; 0.1; 0.2; 0.25; 0.3; 0.35 and 0.5 
GeV. 
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